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Abstract. We present results from simulations to examine how the separation between two con- 
vectively unstable layers affects their interaction. We show that two convectively ustable layers 
remain connected via the overshooting plumes even when they are separated by several pressure 
scale heights. 
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INTRODUCTION 



Some main-sequence stars, such as A-type stars, have a c omplex internal structure in 
whic h there are believed to exist multiple convection zones (|Toomre et all iKupkal 
2005|) . Such a complicated internal structure arises because of the change in the chem- 
ical composition as the distance from the centre of the star is increased. The outer- 
most convection zone, located just below the photophere is belived to be caused by 
the partial ionization of hydrogen and the single ionization of helium. The lower con- 
vection zone is believed to result from the second ionization of helium and is sepa- 
rated from the upper convectio n zone by a re latively narrow convectively stable region 
(|Latour. Toomre & Zahnlll976l : lKupkall2005|) . 

The presence of two convectively unstable zones in such stars raises a number of inter- 
esting questions. However, in this short article we choose to focus on the principal ques- 
tion concerning the degree of intereaction between the convectivley unstable regions. 
Such an issue is motivated by the fact that in stars, without the presence of rigid bound- 
aries, ascending and decending convectively driven motions overshoot the bottom of the 
layer that is convectively unstable. When there are multiple convection zones, in close 
proximity, in stars this overshooting behaviour leads to enhanced transport and commu- 
nication between the unstable layers. Further, due to the close proximity of the unstable 
layers, fascinating dynamics may occur. For example, overshooting plumes form each 
unstable layer could overshoot so far that they continue through the stable region and 
pierce the other convection zone, which would leed to transp ortation of ' contaminants ' 
directly from one convectively unstable region to the other (|Tobias et all ll998L 1200 
In this article we will present some results from simulations to explore how increasing 
the separation affects the connection between the layers. 

It has been conjectured, from analytic theory, that two pressure scale heights 
should be sufficient t o 'di s connect' the two convectively unstable regions 
(ILatour. Toomre & ZahnL 1 19761 : Too mre et al. Here, we will briefly show 

that this is an underestimate and that to 'disconnect' the two convectively unstable 
layers many more scale heights of separation is required. A more detailed exploration 



of this issue can be found in lSilvers & Proctorl (I2007|) . 



MODEL 



Here we examine the evolution of a compressible flui d in a Cartesian layer. The govern- 
ing equations are (IMatthews. Proctor. & Weisslll995h : 
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where z is taken downward, 6 is the dimensionless temperature difference across the 
layer, i?* is the gas constant, m is the polytropic index, K = K / dpQCp^(R*To) is the 
dimensionless thermal diffusivity, y is the ratio of specific heats, T is the stress tensor, 
P = pT and a is the Prandtl number. These equations above are solved using a parallel 
hybrid finite-differenc e/pseudo-spectral code; the most co mprehensive description of the 



code can be found in (Matth ews. Proctor. & Weissl . 11995) 



To obtain a suitable basic state for this problem we allow the thermal profile to be 
non-linear and we take 



K 



K]_ 
2 



K2 +K3 K3 
1 + ^ + ^ tanh 



- f^tanh 



+ 1 ) tanh (— 



(4) 



where A is the characteristic size of the transition region between each of the layers. In 
this work the characteristic sizes of the transition regions are taken to be the same for 
simplicity. The static density and temperature profiles are found by solving the equations 
of hydrostatic balance. To this static state, throughout the domain, random perturbations 
are introduced, with amplitudes which lie within the interval [-0.05,0.05]. 

The aspect ratio for the computational domain in this study is 8:8: where & is 
the total depth of the box, and the domain is assumed to be periodic in x and y. The 
conditions on the upper and lower boundaries are: 
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1 Note that the more detailed exploration of this topic in dSilvers & Proctoit 120071) is in a, smaller, 4:4:5? 
domain. 
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FIGURE 1. The initial temperature and density profiles when 3 = 3.0. 



The system we study has a large number of dimensionless parameters, making it 
impractical to conduct a complete survey. Thus a number are held fixed: o = 1.0, 
mi = = 1.0, = 10.0, 7= 5/3 a nd R a = 1.7xl0 5 . Note that in this paper we use a 
subscript 1 on quantities relating to the upper convection zone. Similarly, subscript 2 will 
be used to denote quantities for the convectively stable layer and 3 to denote quantities 
in the lower convective zone. 

The stiffness parameter, S, provides a useful m easure of the relativ e conductivities 
in this problem (for more detailed discussion see iTobias et al. 52 and S3 are 

related to the various polytropic indexes that appear in the problem via the definition 
that 5 2 = (m 2 - m ad ) /{m ad - mi) and S3 = (m 3 - m ad )(m ad - mi). 



THE EFFECT OF VARYING THE THICKNESS OF THE STABLE 

LAYER 

We begin by considering the case where all three zones have equal depth and so 3 1 = 3.0. 
The static state profiles for this fiducial model are shown in Figure [IJ Note that for this 
fiducial case there are 1.60 pressure scale heights across the convectively stable mid- 
layer, which is rather less than the two pressure scale heights estimate suggested by the 
earlier analytic theory as necessary for true separation. 

The state shown in Figure Q] is perturbed and allowed to evolve to a statistically 
steady state as shown in Figure [2] As one might expect the convection is noticeably 
different in the two convection zones. While Figure [2] provides a useful picture of the 
convective state it is difficult to gauge the motions within the convectively stable region 
from such a figure. To obtain a clearer picture of the potential interactions between the 
two convectively unstable regions it is helpful to note that, if the two layers are to be 
considered as 'independent' from each other, then there needs to be a region between the 
layers where the velocity becomes very small. Therefore, to facilitate a clearer picture 
of the degree of interaction of the layers we calculate the variation in z of the horizontal 
averages of the modulus of the z component of momentum. 

Figure Oa) shows that while there is clearly more vertical motion in the convectively 
unstable regions, as one might anticipate, there is still a non-negligible vertical com- 
ponent of momentum in the middle of the box. Thus the two convection zones are con- 



FIGURE 2. Plot at t = 15.02 Sides of the box show the vertical momentum flux and the top of the box 
shows the vertical momentum flux near the top of the box for the case where 9 = 3.0 and S2 = 5.0. 
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FIGURE 3. Left: Plot at t = 15.32. The variation of the vertical component of velocity as a function 
of height at the centre of the box for the case where & — 3.0 and S2 = 5.0. Right: Plot at t = 15.09. The 
variation of the vertical component of velocity as a function of height at the centre of the box for the case 
where 9 = 5.0 and S 2 = 5.0. 

nected in this case and so there is a conduit for mixing between the two convection zones 
for this level of separation. Any reduction in the width of the convectively stable region 
will decrease the number of pressure scale heights of variation across the layer and in- 
crease the level of interaction between the two unstable zones. Such vigorous motion 
in the convectively stable layer was discussed in the context of the downward directed 
hexagon case of Latour et al. and is further confirmed by a plot of the vertical component 
of velocity as shown in|3^a). 

The most obvious way to limit the interaction of the unstable regions is to increase the 
layer depth, S>, and so increase the width of the stable region if the convectively unstable 
regions have the same height. We have gradually increased the layer depth up to 5.0. This 
corresponds to a pressure scale height variation of 3.77 across the stable layer, which is 
abover the estimate suggested by the earlier analytic theory as the transition between 
connected and unconnected convection layers. 

As for the case when 3> = 3.0, we consider the variation of the planar average vertical 
component of momentum throughout the domain in the established statistically steady 
state, for each of the increased box heights. Figures Ob) show this quantity for each case 
when Ql = 5.0. This figures clearly illustrates that, while there is a significant reduction 



in vertical motion in the convectively stable region, there is a non negligible vertical 
component of momentum throughout the box even for a separation of 3.77 pressure 
scale heights. 



CONCLUSIONS 

In this article we have shown that there is no point at which the vertical component 
of momentum becomes very small even when the number of pressure scale heights of 
separation between the two convectively unstable layers is as large as 3.77. Thus there 
is a route for mixing of passive and dynamic quantities between the two convectively 
unstable regions in all the cases we considered. Only for a substantial increase in the 
box height will the vertical component of momentum fall to a very small value at one 
plane in the box. To find the minimum box height for this to occur would be extremely 
computationally expensive and also uninteresting physically. In A-type main sequence 
stars the separation between the unstable layers does not extend past a couple of scale 
heights. Therefore, on the basis of the present work we can conclude that in A-type stars 
there is a clear connection between the convectively unstable zones that lie immediately 
below the stellar photosphere. 

One of the major objectives of this work is to provide a solid hydrodynamical basis on 
which more complex models can be constructed to understand fully the dynamics that 
occur below the surface of A-type stars. With such a simple model we are not yet in a 
position to address important secondary questions such as the influence of the observed 



chemical anomalies (see discussions in iMichaudl (| 19701) or iVauclair & Vauclairl (|1982t) 
for more details). However, we believe our model provides a platform on which we can 
build, so as to address the effects of rotation, magnetic fields and other issues related to 
the dynamics in these stars and will be the subject of future work. 
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